AN EQUIVALENT OF KRONECKER'S THEOREM FOR 

POWERS OF AN ALGEBRAIC NUMBER AND 
STRUCTURE OF LINEAR RECURRENCES OF FIXED 

LENGTH 

NEVIO DUBBINI AND MAURIZIO MONGE 

Abstract. After defining a notion of e-density, we provide for any in- 
teger m > 1 and real algebraic number a an estimate of the smallest e 
such that the set of vectors of the form it, ta, . . . , ta m ~ 1 ) for t € K is e- 
dense modulo 1 in terms of the multiplicative Mahler measure M{A(x)) 
of the minimal integral polynomial A(x) of a, which is independent of 
m. In particular, we show that if a has degree d it is possible to take 
e = 2Wy/M{A(x)). 

On the other side we show using asymptotic estimates for Toeplitz 
determinants that we cannot have e-density for sufficiently large m when- 
ever e is strictly smaller than l/M(A(x)). In the process of proving this 
we obtain a result of independent interest about the structure of the 
Z-module of integral linear recurrences of fixed length determined by a 
non-monic polynomial. 



1. Introduction 

The classical Kronecker's Theorem in diophantine approximation, in one 
of its different versions, says that if 9i,...,6 m are real numbers linearly 
independent over Q, then (t6i, . . . ,tO m ) is dense modulo 1, i.e. if we call 

n m : M m -> R m /Z m = T m 

the canonical projection, then the set {7i m (t9i, . . . , tO m ), t G M} is dense in 

A question left to investigate by this result is which amount it fails 
for particular choices of which are not linearly independent over Q, 
but have a special form. In this paper we will consider the case where 
the (6i) are powers of an algebraic number a, and more in general when 
the 9i statisfy a linear recurrence relation with characteristic polynomial 
A(x) of degree d, and 9\, 62, ■ ■ ■ , 9d are linearly independent over Q (this 
condition is automatically satisfied when 9i = a* -1 , taking the minimial 
integral polynomial of a as A(x)). 

Before stating our main result, we shall give a definition of e-density and 
recall the notion of Mahler measure. 
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1.1. e-density. For e > and positive integer to, let 

e e" 
~2' 2J 

be the cube with edge lenght e in the torus T m , which we assume centered 
around the origin for convenience. 



C t = n m (I e ) C T m , where I e 



Definition. We will say that a set S C T m is e-dense if S + C ? — T m for 
each e > e, or equivalently if S + C t is dense. A set T e M m will be e-dense 
if T + J f = K m whenever e > e, or equivalently if T + I e is dense. 

Alternatively, it is possible to consider on T m the distance 

d 0O (x,y) = m£^\x-y\ oo , x, y e R m , n m (x) = x, ir m (y) = y|, 

and since the p/2- neighborhood of S with respect to <i 00 (-,-) is precisely 
S + C p for each p > 0, we have that S* will be e-dense if and only if its 
e/2- neighborhood is all of T m , for each e > e (or, equivalently, if its e/2- 
neighborhood is dense). 

1.2. Mahler measure. Suppose given a polynomial of degree d with com- 
plex coefficients 

d d 

M x ) = ^2 aixi = ad ' n ( x ~ ' 

i=0 i=l 

and such that a ^ 0. The Mahler measure of A (re) is defined as 

d 

M(A) = M(A(x)) = \a d \ • JJmaxjl, \oii\}. 

i=i 

The Mahler measure is a notion of complexity that for the minimal 
polynomial of a rational number p/q with (p, q) = 1 reduces to max(|p|, 
Furthermore, each coefficient a« can be written as the sum of (J monomials 
with coefficient 1 in a^, cxi, ■ ■ ■ , a d , each appearing at most once in each 
monomial, and consequently we have that Oj < (^) • M(A), for each i = 
0,...,d. 

Note also that 

i i d 

\ad\ 



M(A(x/2)) = L|t.IJniax{l, 12^-1} 
j'=i 

= |a d | • JJmax |i, |aj|| , 



while 



2- d M(A(2x)) = \a d \-l[m^{l,\^-\} 

3=1 

= \a \ • JJmax jl^T 1 ! , I J = M(i(x/2)), 
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where we called A(x) the conjugate polynomial Yli=o aixd 

We have that M(A(x/2)) > 2-W 2 ^M(A(x)) if A(x) has at least d/2 roots 
> 1 in absolute value, while 2" d M(A(2x)) > 2- [d/2] M(A(x)) if A(x) has at 
least d/2 roots < 1 in absolute value. Since at least one of these conditions 
must hold, we have 

r i 2 d \ 2^ 2 ] 

mm \M(A(x/2))' M{A{2x)) J - M(A(x)) ' 

1.3. Main results. We can now state our main results: 

Theorem 1. Let m > d > be positive integers, and let 0±, . . . , m be real 
numbers such that 0±, . . . ,0^ are linearly independent over Q, and 9d+i, ■ ■ ■ 
are inductively defined by the linear recurrence determined by a primitive 
integral polynomial A(x) of degree d and with nonzero constant coefficient. 
Then 

Se = S e , m = {7r m (*M02, • • .,t6 m ), teR}CT m 
is e-dense for each e at least 

r i 2 d \ 2^ 2 ] 

111111 \M{A{x/2)Y M(A(2x)) J " M{A{x))' 

In other terms, this theorem says that if e is as stated, then for arbitrary 
real numbers Xi, . . . ,x m it is possible to find a real number t and integers 
Pi , . . . , p m such that 

Xi < tOi — pi < Xi + e, for 1 < i < m. 

It is clearly possibile to take as the powers (a 4-1 ) of an algebraic number 
a, and as A(x) as the minimal integral polynomial of a. For each primitive 
integral polynomial A(x) it is also possible to take as (OA any 'sufficiently 
generic' recurrence sequence determined by A(x), such that 9i,...,9 d are 
linearly independent over Q. 

Some computational evidence actually makes us advance the following 

Conjecture 1. In TheoremU\we have e- density for each e at least l/M(A(x)). 

This conjectured optimal density coefficient is not very far from what 
obtained in Theorem dj but it would be optimal because of the following 
negative result: 

Theorem 2. Let 0i,02,- ■ ■ and A(x) be as in Theorem^ Then for each 
e < l/M(A(x)) the set Se, m is not e-dense for m big enough. 

Note that if could be to ruled out the existence of Se )jn which are not 
e-dense for some e < 1 when A(x) is not a product of cyclotomic polynomi- 
als, then a positive solution of Lehmer problem [3 Chap. 1] would follow, 
because the existence of a sequence of polynomials with Mahler measure 
approaching to 1 from above would be automatically ruled out. 

During the proof of this Theorem a result about the structure of the 
module of linear recurrences of fixed length is obtained, it is stated as The- 
orem |3] after some notation has been introduced. 
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2. Preliminaries 

Suppose given a polynomial with complex coefficients 

A{x) = a>dX d + • • • + a±x + a , 

and suppose that ad and ao are both nonzero. For each £ > 1 let's consider 
the rectangular £ x (£ + d) matrix 



def 



(a Q ax ... a d 

a ai ... ad 
do a\ . . . 



ad 



\ ao ai 



\ 



a d J 



lA}e is also the multiplication matrix for sections of power series by the 
polynomial A(x), in the sense that if 

oo oo 
i=0 i=0 

are complex power series such that G(x) = A(x) ■ F(x), then 

/ fd+e-i\ 

/9d+e-i\ 

= {Ah ■ 



9d+l 

V 9d J 



fd+l 
fd 



\ fo I 



and in fact the range of coefficients 0, 1, . . . , d + £ — 1 of F(x) uniquely 
determines the coefficients d, d + 1, ...,d + £— 1 of G(x). 
Moreover, it is easy to see that if A(x) factors as 

s t 
A{x) = B{x) ■ C(x), B(x) = h ix\ C{x) = c i x \ 

then we have 



i=0 



i=0 



IMt = (Bh ■ (Cje+s = \C\t ■ \B\ih 
as it is easy to verify, for instance observing that multiplying a power series 
by A{x) and discarding the lowest d coefficients gives the same result as 
multiplying first by C(x) and discarding the lowest t coefficients, and then 
by B(x) and discarding the lowest s coefficients (and vice- versa). 
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For each polynomial A(x) and positive integer m we will also need the 
lower triangular m x m matrix defined as 

( a d \ 

ad-i ad 



clef 



ao ai . . . ad 

ao ai ... ad 

V a ai 



ad J 



Note that {AJg is embedded in {A}g +d as the last £ rows. Similarly to the 
previous case, we have that if A(x) = B(x)-C(x), then {A} m = {B} m -{C} m . 
In particular, if C(x) splits in linear factors as 



C(x) = ft (x - 7l 



8=1 



we have that 



{C} m = Y\_i X ~ 7i}r 



i=l 



3. Proof of Theorem [Q 
We can now prove Theorem [TJ 

Proof of Theorem [IJ From Kronecker's Approximation Theorem [21 pg. 53- 
54] the set Se, m is dense in the closed subgroup of the torus of the elements 
satisfying the same integral relations as the 9\, . . . , 9 m . 

Let A(x) = Yli=o adxd = ad Yli=i( x ~ a i) be the primitive integral poly- 
nomial defining the recurrence sequence 0\, 6 2l ■ ■ ■ ■ If ko, . . . , fc m -i is an in- 
tegral linear relation satisfied by the 6i (in the sense that fa-ify — 0), 
and K(x) = YlJ^o ^i x% ls ^ ne corresponding polynomial, we can consider 
the remainder R(x) = ^2 i>0 riX l of the division of K(x) by A(x). It is a 
polynomial of degree < d — 1 and its coefficients ro,ri,... still define a 
linear relation of the Q{. But 9\, . . . , 9d are linearly independent over Q, so 
R(x) must be zero. 

Consequently K(x) must be a multiple of A(x) in Q[x], and thus in Z[x] 
by Gauss lemma, being A(x) primitive. This shows that for I = m — d the 
rows of the matrix \A\g generate the Z- module of the integral relations of 
9\, ■ ■ ■ , 9 m . 

The closure of Se, m is thus the projection in T m of the hyperplane 

P = P{A) = jt> G lR m : [A]* •« = ()}, 

the kernel of the linear application \A\i : IR m — > M. e . 

It's inverse image under the projection n m is P + Z m , and since we just 
showed that the sublattice of Z m generated by the row of the matrix {AJg 
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is the saturated sublattice of the integral linear relations of 6\, . . . , 6 m , we 
have that P + Z m is also equal to 

Q = Q(A) = {veR m : {A} e ■ v e Z 1 }. 

In fact, the Z-linear map {AJ, : Z m — > 7/ must be surjective (this can also 
be seen directly considering the t x I minors modulo p for each prime p [H 
Lemma 2, Chap. 1], which have gcd 1 because A(x) is primitive). Now, 
P + Z m C Q clearly, and for each vector w G Q there exist a vector z G Z m 
such that [Afl^ • w = {AJi ■ z, and consequently the difference v = w — z is 
in P, and we have that w = v + z£P + Z m . 

We are now reduced to prove that Q(A) is e-dense in M m , for each e 
satisfying the statement of the theorem. Applying the involution 

, Xi-, ■ ■ ■ , X m — X m ) I ^ (x m , X m —i, . . . , X2, Xl) 

to we obtain a set which clearly has the same e-density properties 

as Q(A), and which is Q(A), where A(x) is the conjugated polynomial 
A(x) = a x d + ■ — h a d _ix + Consequently, we can just prove that Q(A) 
is e-dense for 

e > l/M(A(x/2)) = (\a d \ ■ JJmax |i |^ 

and the e-density for e > 2 d /Af (A(2x)) = l/M(A(x/2)) will follow from the 
same estimate applied to Q(A). 

To obtain this estimate, note first that Q = Q(A) will be e-dense if and 
only if the map 7r^ o {AJe to T l is surjective when restricted to I e = J e>m = 
[— e/2, e/2] m . In fact, Q is e-dense if and only if an arbitrary vector v G IR m 
is contained in Q + I t , and 

(v + I e )nQ ^0 VwGR m 

^ (w + [4-/ e )nzVi VweR' 

applying the matrix {AJt to the expression, and where we denoted {AJi ■ I t 
the image of I € under the map {AJg. This passage must be justified because 
the matrix \A\^ clearly does not have rank m, but since Q contains all 
the vectors that are mapped to Z £ the first intersection will be non-empty 
whenever the second one is (the other direction being trivial). 

We will now factor A[x) as a product of polynomials B(x), C(x) G M.[x], 
with a consequent factorization of [AJf, and prove two different estimates 
for each of the two factors B(x) and C(x). We will select a posteriori the 
factorization which provides the best compound estimate. 

More precisely, let B(x) = Y^ii=o ^ iXi anc ^ C(x) = Ylt=o c i x ii w itli A(x) = 
B(x)-C(x), and recall the induced matrix factorization {A}i = {Bji - [C]f +S . 
To prove that 

n e o lA} t = ix, o lB\ t o {Cj e+S 
is surjective from J e m to T e , we can just prove that 717 o {Bji is surjective 
from Is,e+s for some 6, and that [C]^ +s ■ I ej7n contains Is,i+ S - We anticipate 
that C(x) will be choosed monic and with all roots < 1 in absolute value. 
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3.1. Estimate for B(x). We can take 5 = l/b - Let's show that the 
image the of cube Is,e+ S under \B\g assumes each value modulo Z £ , starting 
with an arbitrary v = (v±, . . . ,vg) G M. and building inductively a vector 
w = (wi, . . . ,wg +s ) in Ig,e+s such that {Bjg ■ w — v G 1} . Suppose that 
w G Is,e+s is suc h that the components with index > i of {Bjg ■ w — v 
are in Z, and observe that while varies in the interval [—5/2,5/2] the 
i-th component of \B\p ■ w varies in an interval large b$5 = 1, while all 
components with index > % of \B\i ■ w stay fixed. Consequently we can 
move Wi in the interval [—5/2, 5/2] to ensure that all components with index 
> i of \B\g ■ w — v are in Z, and since repeating this procedure allows us to 
construct a w with the required properties, our claim follows. 



3.2. Estimate for C(x). IfC(x) = 
that it is possible to take 



i=0 c i x 



[l*=i (a; - 7i) , we will see 



5- 



1 1 



7i 



To prove the assertion, we proceed in the following way. Rather than work- 
ing with [C]^ +s we will work with the nonsingular square matrix {C} m , and 
if we prove that the image under {C} m of I em contains I$ j7n our claim will 
follow, since the image under \C\ m is a just projection on the last t + s 
coordinates of the image under {C} m . 

Now, rather than proving that the image under {C} m of J e ,m contains 
h,mi it w ih be easer to prove that the image under {C}" 1 of Is, m is contained 
in I t m . In particular, {C} m factors as 



{C} m = Y[{X - 7i}m, 



i=i 



and the inverse of a matrix of the form {x — 7} m , for 7 G {7j}i<«<t, is easily 
computed as 



/ 1 



\ 



"7 



"7 



"7 



V 



-7 V 



( 1 




7 


1 


1 


7 






1 


7 2 







\ 



1 

7 



7' 



1 

7 



1 / 

has all 



This shows that if a (possibly complex) vector v 
components with absolute value < p for some real number p > 0, the vector 
w obtained applying the matrix {x — 7}" 1 will have components of the form 



Wr+l 



i=0 
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for some < r < m, and their absolute value can be estimated as 

r 



i=0 



<S|7 i |-kwl<EM , ^iTu/'- 

i=Q i=0 1,1 



Since p is arbitrary we have applying iteratively {x — 7«}m 1 for i = 1, . . . , t 
that the set of complex vectors with all components < 8 in absolute value 
is mapped by {C}" 1 to complex vectors with components of absolute value 
at most e, and consequently Is,m is mapped into J €)Tn , being {C}" 1 a matrix 
with real entries. 

3.3. Conclusion. Let A(x) = B(x)C(x) be a real factoriziation of A(x), 
with 

s s t t 

i=0 i=l i=0 i=l 

and with C(x) monic with all roots < 1 in absolute value. The above 
estimate shows that Q is e-dense for each e at least 

1 3 \ JL i 



and consequently for e at least 



1 fl 11/ 

1^1 ■ II min | R ' t^-\\ = ladl ■ II max {N' 1 - N} 

I I j_Q VII \ 1=1 

since we can take as 7$ precisely the a, t with absolute value < 1/2, and as 
/3i the remaining roots of A[x) (note that C(x) will have real coefficients). 
Since this last expression is clearly not greater than \/M(A(x/2)), the proof 
is complete. □ 



4. Integral linear recurrences of fixed length 

Let A(x) = Y^t=o a i x% ^ e a primitive integral polynomial of degree d with 
ao 7^ 0, and for m > diet A m be the Z-module of the integral vectors in Z m 
that form a recurrence sequence determined by A{x). It is also the module 
of the integral vectors that are in the kernel of the matrix [A] m _^, and as 
such has rank d. In this section we provide some results about the structure 
of this Z-module, and in particular we show that for each prime p we can 
chose a basis with special properties with respect to the p-adic valuation of 
the vectors. 

Only a small amount of these results is required to prove Theorem [21 
but we will state and prove them in full strength nevertheless. We will 
consider the Z p -module Am = A m £g>^ Z p over the ring of p-adic integers 
Z p for simplicity, which is also clearly equal to the vectors in Z™ killed by 
[v4] m _ rf , having [A] m _ d integral coefficients. All considerations can then be 
lifted to Z with arbitrary approximation with respect to the p-adic absolute 
value. 
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Let Vp(-) be the p-adic valuation, and let's recall that the Newton polygon 
of A(x) is the boundary of the polygon in M 2 obtained as the upper convex 
envelop of the points with coordinates 



Suppose that the polygon has r edges with different slopes which connect the 
consecutive pairs of points (w k , v p (a Wk )) for = wq < w\ < ■ ■ ■ < w r = d, let 
— w k — Wk~i be the length of the horizontal projection of the k-th edge, 
and let a k = (v p (a Wk ) — v p (a Wk _ 1 )) / £i be its slope. We have that cr, < Oj 
whenever 1 < i < j < d. 



Newton polygon of A(x) = 9x 4 — 3x 3 — 9x 2 — 2x + 3 for p = 3. 

Let's consider a splitting field K of A(x) over the field of p-adic numbers 
Q p , and let's extend to K the p-adic valuation and absolute value. Call as 
usual Ok the integral closure of Z p in K, which is the local ring formed by 
the element with non-negative valuation, with maximal ideal m generated 
by the uniformizer ir. As it is well known, the polynomial A(x) has precisely 
£k roots with p-adic valuation equal to —<Jk, for 1 < k < r. We will show 
that the slopes of the Newton polygon, and consequently the valuations 
of the roots of A(x), are reflected in the Z p -module Am as stated in the 
following: 

Theorem 3. For each prime p there exist a unique basis of Am such that if 
we form the dx m matrix M = (Mjj) with the basis vectors as row vectors, 
then 

(1) the left dxd submatrix (M it j)i<i<d (resp. the submatrix [M^f] \<i<d 

l<j<d ' m-d+l<j< 

is a block upper (resp. lower) triangular matrix, with blocks Bi, . . . , B r 
(resp. Ci, . . . , C r ) on the diagonal; the number of blocks r is equal 
to the number of edges of the Newton polygon of A{x), and B k and 
Cfc are square matrices of size equal to the lenght £ k of the k-th side 
of the polygon for each 1 < k < r; 

(2) if 1 < s < r is the smallest integer such that o~ s > (or we may 
require alternatively that s is the smallest such that o~ s > 0) then 
Bi, ... , £> s _i, C s , . . . , C r are the identity matrix, and for 1 < k < 
s — 1 the matrix C k has determinant with valuation —a k i k {jn — d), 
and for s < k < r the matrix B k has determinant with valuation 
cr k £ k (m - d); 

(3) walking on the right from the element Mjj = 1 contained in the 
diagonal of the k-th identity block B k for 1 < k < s — 1 (resp. 
walking left from an Mj +m -d,j = 1 contained in the diagonal of the 



(0, v p (a )), (1, Up(ai)), . . . , (d, v p (a d )). 
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identity block Ck, for s < k < r) by t steps we find elements with 
p-adic valuation at least —to~k (resp. at least ta^). 

In other words, the matrix M can be taken with the following form: 



1 






B s 



I3 r 



\ 



B, 



Cn 



\ 



Ci 



Cs-1 






1 



and the behaviour of the p-adic valuation on the rows is controlled by the 
slope of the segments of the Newton polygon of A(x). 

For instance, if A(x) = 9x 4 — 3a; 3 — 9x 2 — 2x + 3 and we take p = 3, then 
for recurrence length m = 10 the matrix M is 



/ 


1 


480 


4203 


3861 


2511 


243 


729 











887 


16853 


33706 


33706 


16853 


33706 







-7722 


-5049 


-3339 


-76419 


33378 


-51543 


1 










887 


16853 


33706 


33706 


16853 


33706 







729 
887 


-44658 
16853 


42822 
16853 


-27306 
16853 


19179 
16853 


3489 
16853 





1 





\o 








729 
1078 


243 
1078 


405 
539 


675 
1078 


423 
539 


48/49 


1 



All denominators are prime with 3, the 3-adic valuations of the elements 
are 



/ 

oo 
oo 
oo 



1 

3 
6 

oo 



2 
3 
3 

oo 



3 
2 
3 
6 



4 
2 
2 

5 



6 
3 
1 

3 



oo 


oo 
2 



oo 
oo 



1 



oo 
oo 
oo 




/ 



and as it is possible observe the valuation increases going right in the first 
row, slowly increases going left in the second and third rows, and increases 
going left on the last row. 

Proof. Let N = (Nij) be the rational dxm matrix obtained putting iVy = 
5{j for 1 < i,j ' < d, and inductively defining the remaining elements in 
each row to form a linear recurrence determined by A(x). We will now 
prove that the square matrix N% = (^Vj,^)i<i,j'<d ^ s nonsingular for £ = 
(1,2, ... ,w,m — d + w + 1, ... ,m), for all w = w , . . . ,w r which are the 
ordinate of a vertex of the Newton polygon of A(x). 

If A(x) has distinct roots ai, . . . ,a d the matrix iV is given by V^ 1 ■ L, 
where 

V = (af^KijXd, L = (c^r 1 ) l<i<d ■ 

l<j<m 

To obtain a formula for the determinant of N% valid for general ccj, let's work 
over C and suppose for a moment that j^A(x) = Yli=i( x ~ a i) where the 
are algebraically independent over C. The determinant of N$ = (-/Vj,^)i<i,j<d 
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is equal to the determinant of = (-^i,^)i<ij<d divided by det V, and this 
turns out to be the Schur function s\ (see [10J) associated to the partition 
A defined as \d-i+i = & — i + 1 for 1 < i < d, evaluated in a%, . . . , ad- 
Note applying the definition of N by linear recurrence that the entries of 
N are polynomial functions in the elementary symmetric functions of the 

and hence they are polynomial in the ctj, and the determinant of each 
submatrix is polynomial in the ai as well. While the «j vary outside of the 
closed algebraic set defined by rL>j( a « — a j) = the determinant of N% is 
equal to the polynomial s\(ai, . . . , ad), and consequently the equality must 
hold for each value of the a,. 

If £ is defined as above, A has precisely d — w parts all equal to m — d, and 
its conjugate partition A' is formed by m — d parts equal to d — w. Recall 
now Jacobi-Trudi's identity 

s\ = o\et(ey_^i + j)i<i t j<k, 

where the are the elementary symmetric functions, and which holds for 
each k at least as big as the number of parts of A'. We are reduced to 
prove that a matrix of the form {ed-w-i+j)i<i,j<m-d is non-singular. After 
evaluation in the roots we have e« = (— l) l dd-i/ad, and flipping the sign of 
rows and columns of even order we can consider the determinant of 

U iflw+i—j)l<.i,j<.m—d 

up to possibly a sign, and discarding a factor a^ m+d . Note that all entries 
on the diagonal are equal to a w . 

If w = wq = (resp. if w = w r = d) then the matrix is lower (resp. 
upper) triangular with a w ^ on the diagonal, and consequently non- 
singular. Suppose w = Wk for some 1 < k < r, let p G K be an element 
with p-adic valuation equal to cr^, and let R be the diagonal matrix with 
1, p, p 2 , . . . , p m ~ d ~ l on the diagonal. The matrix 

±.BT 1 .U.R=( ±*±=*£1 ) 

a "> V a ™ / l<i,j<m-d 

has all entries in Ok, and is upper unitriangular when reduced modulo 
m because it has all 1 on the diagonal, and starting from a w the p-adic 
valuation increases at a rate bigger than cr^. 

Consequently N% is non-singular, and we can consider the matrix Q = 
N^ 1 ■ N. This is already sufficient to show that the matrix M must be 
unique, because suppose that we have another matrix M' with the same 
properties, and suppose that the z-th row differs. If (%)i<j<m is the differ- 
ence of those rows we must have Vj = unless Wk < j < m — d + Wk for some 
k, so taking w = Wk we have that (vj) cannot be a linear combination of the 
rows of Q. But Q has rank d, and the existence of (vj) would imply that 
the module of linear recurrences A m ®^ Q has rank > d, which is absurd. 

Suppose now that w = Wk with 1 < k < s, so that the slope Ok is < 
by definition of s. We will show that all entries in the rows with indices 
w k -i + l,...,Wk — l,Wk are integral, and that the valuation of Qij for 
Wk-i < i < Wk and j > i is at least —<Jk(j — 



12 



N. DUBBINI AND M. MONGE 



In fact, let p have valuation equal to Uk and consider the polynomial 

d — ^ d 

B(x) = ^ P - -x* = b ^ e °kH 



Note that the Newton polygon of B(x) is obtained by the Newton polygon 
of A(x) by subtracting the linear affine function f(x) = a k (x — w) +v p (a w )), 
and this has the effect of moving the k-ih side of the Newton polygon to 
lay on the horizontal axis. In particular, the coefficients of B(x) are such 
that b Wk l and b Wk are ^ (mod m), but bi = (mod m) for i < Wk-i or 
% > w k . 

On the other hand, for each Wk-i < i < Wk the vector {p*~ l Qij)i<j< m 
is a linear recurrence determined by B(x), and we claim that all entries are 
in Ok- In fact, suppose this is not the case, and multiply it by the smallest 
power of the uniformizer n required to make all entries in Ok- Some entry 
will be in Ok \ tn, but the first £k = Wk entries will be in m. When reduced 
modulo m, the subvector (p- ?_l <5jj') u , fe _ 1 <j< m -d+ ?i , fc is a recurrence determined 
by the polynomial 

4 4 

i=0 i=0 

and this recurrence of order £k in Ok/w is supposed to have non-zero entries, 
while the first Ik entries are zero. But this is absurd, and the claim is proved. 
The matrix {p*~ l Qi,j) w k -i<i<w k is invertible, because modulo 

m— <2+«jfc_i<j<m— d+Wk 

m it is obtained as the (m — d)-th power of the companion matrix 

/ -Co/Q fe \ 

1 -ci/c 4 
1 -C2/C4 

V 1 -C4-1/C4 ) 

of C(x), which is invertible. Consequently (Qij) w k _ 1 <i<w k , which 

m— d+Wk-i<j<m— d+w k 

is obtained by it conjugating by the diagonal matrix with diagonal 1, p, . . . , p £k ~ 1 
and multiplicating by a factor p~( m - d \ has the same valuation as p-^m-d) ^ 
i.e. equal to —Oktkipi — d). 

In this way we have built the rows from Wk-i + 1 to Wk of M, and 
proved that Ck = (Qi,j) w k ^<i<w k has determinant with valuation 

m— d+w k _i<j<m—d+w k 

-a k £k{m-d), while B k = (Qi^wk.^ijKwk is the identity. For s < k < r we 
can clearly proceed in a symmetrical way, taking Ck equal to the identity 
and proceeding on the left up to Bk which will have determinant equal to 
odk{™> ~ d). 

The matrix we have built selecting at each step the rows from Wk-\ + 1 
to Wk of Q clearly satisfies all requirements for the matrix M. Furthermore 
all rows are linearly independent, so the module they generate over Z p has 
rank d. 
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To prove that the rows generate all of Am observe that they generate 
a Z p - module contained in Am , and suppose the generated module to be 
strictly contained. By left multiplication by a matrix B with determinant 
G Q p \ Z p we can obtain in any case a basis of Am , but the matrix = 
(Mi£.)i<i t j<d for £ = (1, . . . , w s , m — d + w s + 1, . . . , m) has determinant 
1, and B ■ (and consequently B ■ M) would not have coefficients in Z p , 
which is absurd. □ 

Let's now turn to the Z-module A m again. If E is a Z- module and FC£ 
a submodule with index n, we have that the index of F ®i Z p in E £g>z Z p 
is precisely the biggest power of p dividing n. 

Let m C Z d x Q m ~ d be the Z-module of the linear recurrences deter- 
mined by A(x) such that the first d coordinates are in Z, and similarly let 
9m = O m ®z Zp. The matrix iV of the proof of the Theorem is clearly a 
Z-basis of 9 m . 

For each prime p the matrix M is equal to M^-N where = (Mj^ixij^d 
and £ = (1,2, ... ,d), and as it is possibile to verify immediately has de- 
terminant with the same p-adic valuation as a™~ d , because for s < k < r 
the matrix B/, has determinant with valuation equal to a^i^im — d) = 
{v p {a Wk ) - v p (a 1Bk _ 1 )) (m - d). 

Consequently (9™ : Am) is equal to the biggest power of p dividing 
a™~ d , for each p. We have proved: 

Corollary 1. The module A m has index equal to \ad\ m ~ d in O m . 

5. Proof of Theorem [2] 

Before providing the proof of Theorem 121 let's recall a few facts of linear 
algebra. Let (-, •) be the standard inner product on M. n , and u\,...,Uk G M n , 
for k < n. The Gram matrix of the U{ is the k x k matrix defined as 

G(ui, ...,u k ) = ((ui, Uj))i<ij<k, 

and geometrically its determinant is the square of the volume of the par- 
allelepiped formed by the vectors «j. If the can be completed with 
Uk+i, ...,m„ 6 K" to a basis of M. n , and V\, v 2, ■ ■ ■ , v n G MJ 1 are such that 
they form a pair of biorthonormal bases, then we have 

G(ut, u 2 ,..., u n )G(vi,v 2 , ...,v n ) = I, 

(see [SI §66.2, pag. 66-6]). Since G(u\, . . . , u k ) is a minor of G(ui, u 2 , ■ ■ ■ , u n ) 
and its complementary cofactor is G(vk+i, ■ ■ ■ , v n ), if the u\, . . . , u n form a 
parallelepiped of volume 1 (and consequently G(u\, . . . , u n ) has determinant 
1) we have that 

det G(ui, ...,u k ) = det G(v k+h ...,v n ), 

as it follows from the properties of compound matrices (see [BJ Chap. 1, 
§4, pag. 21, equation (33) in particular]) related to what sometimes is also 
called "Jacobi's Theorem" §14.16]. 
It is now possible to give the 
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Proof of Theorem As in the proof of Theorem [H put £ = m — d, and let 
Q be defined as 



Let oux, . . . , U)d 6 IT 1 be any basis of the (i-dimentional the lattice A m of vec- 
tors in Z m which are killed by \A\i (i.e. form a linear recurrence determined 
by A(x)), and let ei, . . . , e m be the standard basis of M. m . 

Suppose by contradiction e < e < l/M(A(x)) are numbers, independent 
of m, such that Q is e-dense. Since e > e we have n m (Q + [0,e] m ) = 
T m , and since Q is the union of integral translates of the parallelepiped 
formed by Ui, . . . ,Ud the map 7r m : W m — > T m must be surjective on the 
parallelepiped II formed by combinations with coefficients in [0, 1] of the 
vectors ee±, . . . , ee m , Ui, . . . , Ud- 



The map 7r m locally preserves the volume and the image of II is all T m , 
so the volume of II must be > 1. On the other side, the volume of II can 
be computed as the sum of the volumes of the parallelepiped formed by 
all choices of m vectors within eei, . . . , ee m , u>i, . . . , Ud, and note that the 
volume of the parallelepiped formed by, say, ee Sl , . . . , ee s , uJ n , . . . , u Tp with 
p + q = m is not greater than 



being q > m — d and e < 1. On the other side the total number of such 
parallelepipeds is ( m ^ d )- 

The theorem now follows if we prove that we can select a special basis 
oji, . . . , cud of the lattice A to ensure that G(u ri , . . . , u rp ) < K(M)M 2m for 
each subset of the and for each M > M(A), for some constant K(M) 
not dependent of m. In fact, volume of IT must clearly be 



and in particular it — > as m — > oo if we take M(A) < M < 1/e, which we 
can do since we supposed e < 1/M(A). 

To construct the required basis, let d, . . . , Q £ Q m be such that the 
j-th coordinate of Q is 5ij for 1 < i,j < d, and defining the remaining 
coordinates by the linear recurrence determined by A(x). The Q are clearly 
a basis of the Z-module O m of vectors in Z d x Qf which are killed by {Ajg, 
and we have by Corollary [1] that (A m : 6 m ) = \a>d\ e . Consequently, a basis 
(cUi) of A m can be obtained by the (Q) applying a matrix W = (Wi t j)i<ij<d 
with determinant \a.d\ e - Changing it by left multiplication with an element 
of SL{d,1i) the matrix W can be supposed to be in Hermite Normal Form 
[51 §23.2, pag. 23-6,23-7], and in particular upper triangular and such that 
\Wij\ < \Wjj\ whenever i < j. 
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We will now show that we can bound the Gramian of a subset of the Ui in 
terms of the Gramians of all subsets of the Q and the determinant of (Wij) 
squared (which we know to be equal to |a<i| 2 ^). In fact, let co n , . . . ,co rp , for 
1 < p < d and 1 < r% < ■ ■ ■ < r p < d, be a subset of the Wj. The quantity 
yjG(bj ri , • • • , w rp ), being the volume of the parallelepiped formed by the co ri , 
is also equal to 

sup (j)(u n ,...,u rp ), 
<t>e<$> P 

the sup being taking withing all the elements of the exterior power A*(M m ) 
of the form 



3>p = { n i A n 2 A ■ ■ ■ A n p , for orthonormal ni, . . . , n p G R m } C A* 



where R m is identified with the dual through (•,•). This can now be esti- 
mated as 

d d 

sup 4>{J2 W ruj (j, W»wCi) 

= sup Y <f>(W ri , sl ( Sl , . . . ,W rp , Sp ( Sp ) 
^^i< Sl ,... lSp <d 

< SU P W ru'l ■■■ W r P ,s P ■<f>((si,---,(s p ) 
l< Sl ,... )Sp <d^ e *f 

< \ w w- w wp\ ■ su p <f>(tsi> • • • Xs P ) 

l<sx,...,s p <d * g *p 

< det(Wy) • SU P ^(C-n • • ■ > Csp), 

since we can discard the summands where Sj = Sj for some i ^ j, and (Wij) 
is upper triangular and with integral entries, 



< dettwy - <F • max JG(C S1 , . . . , &„)■ 

l<si<— <s p <<2 V * 

Note in particular that, for each M > M(A(a?)), this allows us to convert 
an estimate of the form 

G(Cn: ■ • • , Cr„) < K(M) [\a^\ ■ M] 2m , for all 1 < n < - • • < r p < d 

(1) 

into one of the form 

G(u n , oo rp ) < K'(M)M 2m , for all 1 < n < ■ • ■ < r p < d 

for some other constant K'(M) = d p K(M) (thus uniformly in m). 

To prove the ([TJ, let -B(x) = ^(x) (so that M(5) = (a" 1 1 • M(A)), 
and let (-Bj)i<«<£ be the row vectors of the rational matrix \E\^. Note that it 
can be completed to a matrix with determinant 1 inserting the row vectors 
ei,...,ed of the standard basis (ej)i<j< m of R m , and the Q are dual to 
the Ci in the basis ei, . . . , e^, -B 1; . . . , with respect to the standard scalar 
product (•, •). 
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Consequently, let ( ri , . . . , ( rp for 1 < r\ < ■ ■ ■ < r p < d be a sub- 
set of the Q, and let 1 < si < ■ • ■ < s q < d be the complementary set 
of the rj in 1, 2, . . . , d. As observed before the proof, the Gram determi- 
nant of ( ri , . . . , ( rp is the same as the Gram determinant of the vectors 
e Sl , . . . , e Sq , Bi, . . . , Bf . Note that this Gram matrix is a matrix with a very 
special form, and after a complicated part we have that if j, k are > q 
the (j + l,k + 1) entry is equal to the (j, k) entry being (Bj_ g , B k _ q ) = 
(Bj + i_ q , B k+1 _ q ), and in particular it is a perturbed Toeplitz matrix. 

The asymptotic behaviour of matrices that are perturbed Toeplitz ma- 
trices have been deeply studied by Bump-Diaconis [2], Tracy- Widom [12] 
and Lyons [9]. In particular, the following theorem of Lyons [9j Theorem 
3.1 ] suits perfectly our needs. 

Let A the Lebesgue measure on the unit circle, and let E n (t) = e 2nmt . 
Let Poly n be the linear span of E ,Ei, . . . ,E n . Let H 2 (T) be the Hardy 
space of the functions / G L 2 (\) such that all negative Fourier coefficients 
vanish. Let P^i be the orthogonal projection L P (X) — > H 2 (T). For / > 
with log / G L^A), let 



= exp ■ 



1 f E 1 (t) + z 



\ogf(t)d\(t), 



2 J T E^t) - z 
for \z\ < 1. Let <fif{t) the be the outer function 

<P f (t) := lim ^(ri^)), 

which exists A-almost everywhere [TT], and satisfies \(fif\ 2 = f, A-almost 
everywhere. If fi is a measure on T, let its Fourier coefficients be fi(n) = 



It 
j)) 



-2nint 



d/i(t), and let D n be the determinant of the Toeplitz matrix (fl(k- 
o<j,k<n- Then we have the following result: 



Theorem 4 (R. Lyons). Let fi be a finite positive measure on T with infinite 
support. Let f := [dfi/dX] be the Radon-Nikodym derivative of the absolutely 
continous part of ji. 

Given any functions f , . . . , f r , g , . . . , g r G L 2 {^), let Fj := P H ^(fj(f>f) 
and Gj := P H 2(gj(f)f) . Defi 



ne 



and 



We have 



Pi 



Qj 



Ej if - 1 > j > 



9j ifO<j<r, 

Ej if - 1 > j > -n. 



lim D n -i(n) dct 



PjQk dfi 



det 



—n<j,k<r 



Note that we can take n 



q-h n = for / 



FjGk dX 

- 0<j,k<r 

B(e 2vit )B(e- 2mt ) 



and fj = gj = e -^j+^B(e- 27Tit )-\ which clearly G L 2 (/i), for < j < q-1. 
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In this way after a trivial computation we get 

PjQkdfi = G{e si , . . .,e 8q ,B u . . .,B t ), 

J —n<j,k<r 

because for instance for 1 < j < q and 1 < k < £ we have 

Pj-iQljidfi = [ e -^^B(e- 2mt y l ■ e 2wikt ■ B(e 27rit )B(e- 2nit )dt 
Jo 

= [ e 27Ti{ - s J +k)t ■ B(e 2nit )dt 
Jo 

= b S j-k = B k ). 

On the other hand (ft(k — i))o<j,k<n after is the Toeplitz matrix associ- 
ated to the Laurent polynomial B(x)B(x~ 1 ) as can be readily verified, and 
consequently _D„_i(/i) = G(B 1: . . . , B e ). 

The theorem tells us that for £ — > oo the asymptotic value of the ratio 

det G(e Sl , . ..,e Sq ,B u ..., B t )f det G{B X , ...,B t ) 

is finite for each Si, . . . , s g e {1, . . . , d}, so we are reduced to show that 

det G(B X , ...,B e )< K(M)M 2m 

for each M > M(B) for some constant K(M). 

But a nice result on the determinant of banded Hermitiatian Toeplitz 
matrices is available, thanks to the Trench's Formula [H Theorem 2.10, 
pag. 41]: let C(x) = ^Y^j=-r c i x ^ ^ e a Laurent polynomial, and let g n (z) be 
the row 

g n (z) = (1, z, z 2 , . . . , z r ~\ z n+r , z n+r+ \ . . . , z n+r+s - x ). 

Let £i, ■ ■ ■ ,£k be the distinct roots of C(x), and let . . . , fik be their mul- 
tiplicities. We define G n as the determinant of the (r + s) x (r + s) matrix 
r„ whose first /ii rows are g n (£i), g' n (£i), • • • , g^ 1 whose next /i 2 rows 

are g n (£ 2 ), #n(&)> • • • i9n _1) (6), and so on. 

Theorem 5 (Trench). We have Go ^ 0, and 

Dn-tiC) = (-irc^, for every n > 1. 

Applying the Theorem with C(x) = B(x)B(x' 1 ) and n = £, we just have 
to show that c n s G n is smaller than K(M)M n for each M > M(C) = M(B) 2 , 
for some K(M). But if we expand the determinant of the (r + s) x (r + s) 
matrix T n we have that it can be written as sum of monomials in the £j with 
polynomials in n as coefficients, where in each monomial £j appears with 
exponent smaller than (n+r+s)/ij. Consequently estimating each monomial 
the determinant is < P(n)M(j-C) n+r+s for some polynomial P(n) in n, and 
in this way we obtain an estimate of the form c™G n < K(M)M n for each 
M > M(C), as required. □ 
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